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ABSTRACT 

Turbulent  temperature  fluctuations  in  the  atmospheric 
boundary  layer  measured  at  2m,  7m  and  23m  in  Ris0,  Denmark, 
were  analyzed  with  particular  emphasis  placed  on  determining 
characteristics  of  the  high  frequency  region  of  the  spectra 
of  these  fluctuations.   The  shape  of  the  high  wave  number 
one-dimensional  temperature  spectrum  and  an  estimate  of  the 
Kolmogorov  scalar  constant  were  determined.   Comparisons 
of  high  frequency  spectral  regions  of  temperature  and 
velocity  fluctuations  were  made. 
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I.   INTRODUCTION 

Several  attempts  over  the  past  few  years  have  been  made 
to  measure  high  frequency  temperature  variations  in  the 
atmospheric  boundary  layer.   Sensors  with  a  sufficiently 
high  frequency  response  and  spatial  resolution  have  been 
devised  to  overcome  the  inherent  problems  previously 
encountered  in  trying  to  measure  small  scale,  high  frequency 
temperature  fluctuations.   In  addition  to  the  sensitive 
sensors,  low-noise  electronic  equipment  has  become  available 
to  record  and  analyze  small  scale  fluctuations. 

The  object  of  this  research  was  to  utilize  an  analog 
analysis  scheme  to  determine  the  shape  of  the  high  wave 
number  (small  scale)  one-dimensional  temperature  spectrum 
and  to  evaluate  the  scalar  constant,  Kfi ' ,  which  appears  in 
the  expression  for  the  -5/3  form  of  the  temperature  spectrum 
In  addition,  simultaneous  analyses  of  velocity  fluctuations 
were  conducted  in  order  to  obtain  the  necessary  information 
ultimately  required  to  evaluate  Kfi ' . 

A  knowledge  of  the  small  scale  temperature  distribution 

in  a  turbulent  field  has  a  number  of  applications.   It  is 

necessary  when  considering  the  heat  budget  which  is  of 

utmost  importance  in  the  fields  of  oceanography  and 

meteorology.   Sensible  heat  flux  can  be  determined  by 

the  covariance  of  measured  temperature  fluctuations  and 

vertical  components  of  the  simultaneous  wind  velocity 

fluctuations . 
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Temperature  inhomogeneities  in  the  atmosphere  are 
one  of  the  quantities  associated  with  the  scattering  of 
acoustic  waves  and  electromagnetic  radiation.   This 
scattering,  which  is  a  major  concern  in  the  fields  of 
acoustics  and  optics,  is  often  related  to  variations  in 
the  refractive  index. 

Further  applications  of  a  knowledge  of  the  small  scale 
temperature  distribution  are  found  in  the  areas  of  physi- 
cal chemistry  and  thermodynamics. 
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II.   THEORETICAL  DEVELOPMENT 

A.   VELOCITY  SPECTRUM 

Small  scale  turbulence  in  a  fluid  is  a  feature  of 
interest  to  the  environmental  scientist  due  to  its  effect 
on  heat,  momentum  and  energy  transport.   The  one-dimensional 
velocity  (kinetic  energy)  spectrum  is  given  by 

u2  =  S°   <Kk)dk  (1) 

0 

where  u   is  the  total  variance  of  the  turbulent  velocity 
fluctuations  and  k  is  the  radian  wave  number.   The  rate  of 
energy  dissipation  in  an  isotropic,  turbulent  flow  field  is 
given  by 

e  =  15v   /"  k2<Kk)dk  (2) 

0 

2 
where  the  function  k  $(k),  the  dissipation  spectrum, 

statistically  describes  the  rate  of  decay  of  turbulent 

2     - 1 
kinetic  energy  to  heat,   v (cm  .sec   )  is  the  kinematic 

viscosity  of  the  fluid. 

Modern  turbulence  research  had  as  its  impetus  the 

concept  of  a  continuous  cascade  of  kinetic  energy  from 

large  scale  anistropic  velocity  fluctuations  to  smaller 

scales  where  the  motion  became  increasingly  isotropic  and 

at  which  the  energy  is  dissipated  as  heat.   Small  scale 


12 


motions  which  are  characterized  by  small  time  scales  are 
statistically  independent  of  the  relatively  slow  large 
scale  turbulence  derived  from  the  mean  flow. 

This  concept  was  presented  by  the  Russian  mathematician, 
A.  N.  Kolmogorov  in  1941.   Kolmogorov  proposed  the  idea  of 
local  isotropy,  i.e.,  at  some  point  in  the  cascade  of 
kinetic  energy  from  large  scales  to  small  scales,  the 
turbulent  fluid  motion  becomes  isotropic.   This  idea  was 
described  by  Kolmogorov  in  two  hypotheses: 

(1)  If  the  Reynolds  Number  (ratio  of  inertial  to 
viscous  forces)  is  sufficiently  large,  then  there  exists 
a  range  of  small  scales  in  which  statistically,  the  turbu- 
lence is  basically  steady-state  and  is  dependent  only  on  the 

2    - 1 
rate  of  dissipation  of  kinetic  energy,  e[(cm/sec)  .sec   )], 

the  wave  number,  k(cm  )  and  the  kinematic  viscosity, 

2     -1 
v(cm  .sec   ).   Accordingly,  dimensional  analysis  yields 

*(k)  =  (ev5)1/4  F(k/ks)  (3) 

for  the  form  of  the  one -dimensional  energy  spectrum  and 

k2<Kk)  =  ks2(ev5)1/4  (k/ks)2  F(k/ks)  (4) 

for  the  form  of  the  energy  dissipation  spectrum.   Here, 
k  =  (e/v  )     is  the  Kolmogorov  wave  number,  1/k   is  known 
as  the  Kolmogorov  microscale.   F(k/k  )  is  a  universal  func- 
tion which  is  valid  for  all  fields  of  turbulence  provided 
the  Reynolds  Number  is  sufficiently  large. 
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(2)  If  the  Reynolds  Number  is  very  large  (on  the  order 
of  10   and  greater)  then  a  range  of  scales  exists  wherein 
the  energy  cascades  independent  of  viscosity.   In  this 
range  of  wave  numbers,  the  decay  of  the  flow  energy  is 
characterized  by  a  -5/3  dependence  upon  wave  number.   This 
range  in  which  the  energy  is  transferred  inertially  is 
known  as  the  inertial  subrange  and  characterizes  the  turbu- 
lence scales  too  large  to  be  influenced  by  viscosity  but 
small  enough  to  be  separate  from  the  mean  flow.   In  the 
inertial  subrange, 


$(k)  =  K»e2/3  k"5/3  (5) 


where  $ (k)  is  the  one-dimensional  energy  spectrum 

2 
[(cm/sec)  .sec]  and  K1  is  a  universal  constant.   Previous 

researchers  have  reported  values  of  K'  ranging  from  0.48 

to  0.69  although  the  lower  values  are  considered  to  be 

more  accurate. 

B.   TEMPERATURE  SPECTRUM 

It  is  possible  to  apply  the  Kolmogorov  ideas  to 
turbulent  fluctuations  of  passive  scalars  such  as  tempera 
ture .   Temperature  is  a  passive  scalar  when  it  is  trans- 
ported in  a  fluid  but  does  not  introduce  buoyancy  effects 
If  0(C°)  is  the  deviation  of  the  temperature  from  its 
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mean,  the  one  dimensional  temperature  fluctuation  spectrum 
is  defined  as 

e2  =  /°°  *  oo dk  (c0)2  .  (6) 

0    e 

The  rate  of  temperature  dissipation  in  an  isotropic  flow 
field  is 


eQ  =  6k  f"   k2<*>  (k)dk  (C°)2.sec"1  .  (7) 

9      Q      6 

2 
k  is  the  molecular  diffusivity  and  k  $Q (kj  is  the  scalar 

dissipation  spectrum  which  describes  the  distribution  with 

2 
wave  number  of  the  rate  of  decay  of  the  quantity  0  ,  the 

temperature  variance. 

Again,  according  to  Kolmogorov: 

(1)  If  the  Reynolds  number  is  sufficiently  large  there 
exists  a  range  of  small  scale  temperature  fluctuations  which 
statistically  are  dependent  only  on  e,  v,  efi  (the  tempera- 
ture dissipation  rate)  and  k.   Once  again  dimensional 
analysis  yields 

<fre(k)  =  eee"3/4  v5/4  H(a,k/ks)  (8) 

for  the  one-dimensional  temperature  spectrum  and 

k2*e(k)  =  ks2eee"3/4  v5/4  (k/ks)2  H(o,k/k J        (9) 
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for  the  temperature  dissipation  spectrum,  o   is  the 
Prandtl  number  (a  =  v/k)  ,  the  ratio  of  kinematic  viscosity 
to  thermal  diffusivity.   In  a  gas  o      is  approximately  1. 
For  air  at  20°C  the  value  is  near  0.7. 

C2)  If  the  Reynolds  Number  is  even  larger,  then  the 
temperature  spectrum  is  independent  of  v  and  k  and  is 
determined  only  by  e  and  efi.   Dimensional  analysis 
[Corrsin,  1951]  yields 


*fl(k)  =  K  'efie'1/3  k~5/3.  (10) 


K  '  is  supposedly  a  universal  constant  although  previous 
researchers  have  reported  values  ranging  from  0.2  to  2.0. 

In  view  of  this  range  of  reported  values  of  Kfl '  it 
was  the  intent  of  this  research  to  determine  the  shape  of 
one-dimensional  temperature  spectra  at  high  wave  numbers 
and  to  evaluate  K  '  from  the  data  analyzed.   In  this  way, 
a  larger  data  base  could  be  established  from  which  deter- 
minations for  KQ '  could  be  resolved. 
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III.   INSTRUMENTATION  AND  EXPERIMENTAL  APPROACH 

A.  BACKGROUND 

Sensors  and  electronic  systems  with  the  necessary 
sensitivity  and  response  to  measure  high  frequency,  small 
scale  atmospheric  temperature  fluctuations  were  developed 
in  the  early  1970's  by  National  Electrolab  Associates 
Limited  of  Vancouver,  British  Columbia,  in  consultation 
with  the  Institute  of  Oceanography,  University  of  British 
Columbia.   Previous  electronics  allowed  the  temperature 
spectrum  to  be  measured  to  just  beyond  the  inertial  sub- 
range but  not  to  sufficiently  high  wave  numbers  to 
unequivocally  evaluate  the  area  under  the  scalar  dissipa- 
tion spectrum. 

B.  TEMPERATURE  MEASURING  SYSTEM 
1 .   Sensor 

The  sensor  of  temperature  fluctuations  utilized 
was  a  resistance  thermometer  consisting  of  a  platinum  wire 
0.25  ym  (0.00001  in.)  in  diameter  and  0.30  mm  in  length. 
The  wire  has  a  time  constant  of  less  than  10  ys  in  a  4  m/sec 
air  flow.   With  a  detection  current  of  50  yamps,  its 
sensitivity  is  0.063  mv/°C.   The  sensor  is  made  from 
Wollaston  wire  consisting  of  a  silver  jacket  about  a 
platinum  core.   The  outside  diameter  of  the  jacket  is 
45  ym  (0.0018  in.).   The  jacket  and  core  are  prestressed 
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into  a  V-shape  and  then  soldered  onto  the  end  of  a  standard 
hot  wire  anemometer  probe  (Figure  1) .   The  sensor  has  a 
frequency  response  which  is  essentially  flat  to  5KHz. 
2.   Bridge 

Utilizing  a  thin  platinum  wire  to  measure  temperature 
fluctuations  takes  advantage  of  the  fact  that  a  slightly 
heated  wire  is  basically  a  simple  resistance  thermometer. 
In  order  to  sense  a  change  in  resistance  in  the  wire  due 
to  variations  in  the  ambient  temperature,  a  small  heating 
current  is  passed  through  the  wire  and  the  resulting  voltage 
drop  can  be  measured.   If  the  heating  is  kept  low,  the 
resistance  changes  due  to  cooling  by  wind  are  negligible 
compared  to  resistance  changes  caused  by  temperature  fluc- 
tuations.  A  bridge  circuit  is  utilized  to  ensure  optimum 
current  flow  through  the  platinum  core. 

The  bridge  circuit  is  a  modified  Wheatstone  bridge 
consisting  of  two  2.21Kft  resistors,  sensor,  probe  and  cable, 
balance  adjustment  resistors  and  capacitors.   The  fixed 
and  variable  capacitors  (Figure  2)  paralleling  the  balance 
resistance  are  necessary  to  match  the  probe  and  cable 
capacitance.   The  bridge  amplifier  is  a  differential 
amplifier  with  a  fixed  gain  of  60dB  and  a  frequency 
response  flat  to  150KHz. 
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Figure   1.  Wollaston    Wire   Mounted  on   Probe 
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C.      EXPERIMENTAL    PROCEDURE 

1 .  Location 

The  data  were  collected  at  the  meteorological 
field  site  operated  by  the  Meteorological  Group  of  the 
Danish  Atomic  Energy  Commission,  Research  Establishment, 
Ris$.   The  site  features  a  130M  meteorological  tower 
instrumented  at  seven  levels.   In  addition,  miscellaneous 
instrumentation  can  be  mounted  at  arbitrary  levels  provided 
only  that  the  electronics  is  sufficient  to  drive  the  cables 
from  the  instrumentation  to  the  junction  boxes  at  the 
seven  standard  levels.   The  tower  stands  on  the  Ris0 
Peninsula  0.5  Km  from  Roskilde  Fjord.   Signals  from  tower 
instrumentation  are  carried  by  cable  to  a  meteorological 
station  50m  downwind  of  the  tower  where  both  analog  and 
digital  recording  equipment  are  housed.   During  the  course 
of  the  experiments  winds  blew  essentially  from  the  west 
off  the  fjord.   Data  were  not  collected  unless  the  wind 
was  from  off  the  water. 

The  Ris^  Peninsula  is  6Km  north  of  the  ancient 
Danish  town  of  Roskilde  which  in  turn  is  30Km  approximately 
southeast  of  Copenhagen  on  the  island  of  Sjaelland. 

2 .  Equipment 

a.   Meteorological  Tower  Arrangement 

Mean  wind  speed  and  direction,  temperature 
and  humidity  were  measured  at  7,  23,  39,  56,  72  and  96m 
(Figure  3) .   Sonic  anemometer  measurements  were  made  at 
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Figure  3.   Ris^  Research  Establishment 
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2,  23  and  72m.   Fine  structure  turbulence  measurements 
were  made  at  different  heights  during  the  course  of  the 
experiment.   In  general,  Naval  Postgraduate  School  (NPS) 
equipment  was  deployed  at  2m,  7m  or  23m.   Oregon  State 
University  (OSU)  equipment  was  always  at  23m.   University 
of  California,  San  Diego  (UCSD)  equipment  was  placed  very 
near  the  top  of  the  mast  at  123m.   Danish  Atomic  Energy 
Commission  (AEK)  equipment  was  located  at  2m.   Each  group 
deployed  sensors  for  measuring  turbulent  temperature  and 
velocity  fluctuations.   The  main  purpose  of  the  experiment 
was  to  make  simultaneous  measurements  of  atmospheric 
temperature  fluctuations  by  different  systems,  at  the  same 
level  and  at  different  levels.   The  temperature  sensor  in 
each  case  was  a  platinum  wire.   The  diameter  of  the  NPS 
wire  was  0.25  urn,  of  OSU  0.625  ym  and  of  UCSD  both  0.25 
and  0.625  urn.   By  making  these  simultaneous  measurements 
with  similar  sensors  but  different  electronics,  it  was 
hoped  (1)  to  find  in  what  way  they  differed  and  (2)  more 
importantly,  to  resolve  some  of  the  discrepancies  surround 
ing  the  value  of  the  scalar  constant  K  ' . 

The  main  comparisons  made  in  this  thesis  are 
between  the  NPS  and  OSU  systems. 

b.   NPS  Platinum  Resistance  Thermometer 

At  the  2m  height  the  sensor  was  mounted  in  an 
aluminum  holder  which  could  hold  as  many  as  four  such 
sensors.   For  the  data  discussed  in  this  thesis  the 
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temperature  sensor  was  placed  5mm  below  a  Disa  hot  wire 
anemometer.   The  other  two  positions  were  left  vacant. 

At  7m  the  same  aluminum  holder  was  used  but  a 
Thermosys terns  probe  was  used  for  the  hot  wire.   This  probe 
was  vertically  aligned  and  then  rotated  to  be  within  5mm 
to  the  side  of  the  temperature  sensor. 

A  somewhat  similar  arrangement  was  made  at 
23m  except  it  was  more  difficult  to  get  the  temperature 
probe  close  to  the  hot  wire  probe.   There  was  a  separation 
of  approximately  1cm. 

c.   Hot  Wire  Anemometer 

A  variety  of  hot  wire  anemometers  was  used. 
They  were  Disa  (AEK,  UCSD) ,  Thermosys terns  (NPS,  OSU) , 
and  a  new  wind-vane  probe  system  (AEK)  designed  by  Larsen 
and  Busch  (1974)  . 

Disa  probes  were  mounted  in  the  aluminum  holder 
previously  described  whereas  the  other  systems  had  indi- 
vidual clamps. 

3 .   Recording  Systems 

NPS,  OSU  and  UCSD  used  analog  tape  recorders. 
AEK  recorded  digitally  on  magnetic  tape.   The  NPS  data 
analyzed  for  this  thesis  were  recorded  on  an  Ampex  FR  1300 
tape  recorder  at  7  1/2  ips  using  FM  electronics.   At  this 
speed  the  3dB  point  is  at  2.5KHz.   Each  NPS  temperature 
signal  was  recorded  both  direct  and  differentiated  (with 
respect  to  time) .   The  object  of  the  differentiation  is 
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to  increase  the  level  of  the  signal  at  high  frequencies 
above  the  noise  level.   The  differentiator  circuit  had  a 
gain  of  6dB/octave  over  the  differentiating  band  of 
frequencies.   This  band  could  be  varied  from  DC  to  500Hz, 
1000Hz,  1500Hz,  2000Hz  or  3000Hz.   The  frequency  at  which 
unity  gain  occurred  could  also  be  varied.   Beyond  the  band 
of  differentiation,  the  circuit  rolled  off  at  6dB/octave. 
No  other  filtering  of  the  signals  was  done.   NPS  hot  wire 
signals  were  recorded  directly  without  filtering  or 
differentiation. 

The  input  level  of  the  tape  recorder  was  set  to 
handle  1  volt  rms  signals. 
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IV.   ANALYSIS  PROCEDURES 

A.   SELECTION  OF  DATA 

Generally  when  concerned  with  the  analysis  of  high 
frequency  turbulent  temperature  and  velocity  signals, 
it  is  not  necessary  to  examine  long  sections  of  data. 
It  was  determined,  however,  in  the  preliminary  analysis 
of  this  data  that  in  order  to  obtain  characteristic 
spectra,  it  would  be  necessary  to  extend  the  record  length 
to  10  minutes.   This  was  due  to  varying  levels  of  the 
recorded  signals,  necessitating  the  utilization  of  a 
broadband  amplifier  to  maintain  the  signal  level  at  an 
almost  constant  level. 

The  primary  considerations  for  selecting  a  section  of 
data  for  analysis  were  first  that  the  temperature  and 
velocity  fluctuations  be  considered  typical  of  the  record 
and  secondly,  if  possible,  correspond  to  sections  that  were 
being  analyzed  at  Oregon  State  University  for  future  com- 
parison of  results  (Figure  4) .   By  trying  to  adhere  to 
these  considerations,  the  following  data  sections  were 
analyzed:   a  10  minute  length  of  record  from  Run  5  commenc- 
ing at  1537  local  Ris$  time  on  30  August  1974,  and  a  10 
minute  length  of  record  from  Run  6  commencing  at  1657  on 
30  August  1974  (Table  I) . 
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Figure  4.   Example  of  Data  Selected  for  Analysis 
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B.   ANALOG  EQUIPMENT  AND  PROCEDURES 

The  analog  analysis  scheme  consisted  of  utilizing  the 
following  state-of-the-art  equipment  to  obtain  spectra  of 
the  selected  sections  (Figure  5) : 

A  Honeywell  Model  Ninety-Six  Magnetic  Tape  Recorder/ 
Reproducer  System  was  used  to  playback  the  recorded 
signals.   This  device  is  capable  of  recording/reproducing 
any  combination  of  multi-channel  direct  (analog)  and  FM 
data  at  nine  selectable  servo-controlled  tape  speeds 
from  15/16  through  240  inches  per  second  (ips) .   Its 
standard  magnetic  assemblies  are  IRIG  compatible.   It  is 
characterized  by  a  tape  speed  accuracy  of  0.11  and  a 
bandwidth  of  0-2.5KHz  (within  ldB)  at  a  signal-to-noise 
ratio  of  50dB  when  operated  at  7-1/2  ips. 

The  playback  signals  from  the  reproducing  system  were 
fed  to  a  Preston  8300  XWB  Amplifier  with  a  bandwidth  of 
lOOKHz  and  selectable  gains  of  1,  2,  5,  10,  20  and  50. 

The  output  from  the  amplifier  was  fed  into  the  input 
of  a  Federal  Scientific  UA-500  Ubiquitous  Spectrum  Analyzer 
The  selected  analysis  range  was  0-2KHz  at  a  sample  rate 
of  6000hz.   The  analyzer  has  a  3dB  bandwidth  of  6hz  and 
a  frequency  accuracy  of  ±0.21  of  the  full  analysis  range. 

The  input  signal  to  the  analyzer  was  monitored  by  a 
Hewlett-Packard  3400A  RMS  Responding  Voltmeter  and  a 
Tecktronix  Type  502A  Dual-Beam  Oscilloscope. 
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The  spectral  output  was  visually  monitored  on  a 
HP  1300A  X-Y  Display  and  plotted  on  a  Mosely  7035B  X-Y 
Solid  State  Recorder  produced  by  Hewlett-Packard. 

C.   DATA  INTERPRETATION 

Theoretical  work  almost  always  expresses  turbulence 
results  in  terms  of  wave  number  while  actual  turbulence 
measurements  are  spectrally  computed  as  functions  of 
frequency.   These  two  approaches  are  related  by  Taylor's 
hypothesis 


k   =   -y-  (11) 


where  k  is  radian  wave  number  (cm  ) ,  f  is  frequency  (Hz) 
and  U  is  mean  wind  speed  (cm/sec) .   Equation  11  implies 
that  turbulence  with  scales  near  1/k  are  "frozen"  during 
the  period  these  scales  are  swept  past  the  sensor.   The 
maximum  frequency  analyzed  was  2KHz  at  a  mean  wind  speed 
of  5  m/sec  which  indicated  that  a  maximum  wave  number  of 
25cm   would  have  to  be  resolved.   This  requires  that  the 
sensor  does  not  exceed  0.4mm  in  length  and  have  a  fre- 
quency response  essentially  flat  to  2KHz.   Both  of  these 
requirements  were  met. 

Since  Eq.  11  is  linear,  the  shape  of  a  wave  number 
spectrum  is  the  same  as  that  of  a  frequency  spectrum. 
Frequency  spectra  and  wave  number  spectra  are  related  by 
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f<fr(f)  =  k<j)(k)  .  (12) 

The  inertial  subrange  expression  (eq.  5)  in  terms  of  f 
becomes 


<Kf)  =  (g1)    K'e2/3  f"5/3  .  (13) 


The  rate  of  kinetic  energy  dissipation  expression  (eq.  2) 
becomes 


2 
e  =  15v  (,^L)    /"  f2cj)(f)df  .  (14) 

u    0 


The  resulting  expressions  for  the . temperature  spectrum  in 
terms  of  f  are 


2/3 


♦  ft(£)  =  (^)     K  'eQe"1/3  f'5/3  (15) 


and 


2   w 

e0  =  6k  ^    /rof24»Q(f)df  .  (16) 


Equations  13,  14,  15  and  16  were  the  expressions  used 
to  evaluate  e,  eQ  and  KQ ' . 
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V.   RESULTS 

A.   VELOCITY  SPECTRUM  RESULTS 
1 .   Calculation  of  e 

a.   Indirect  Method 

An  indirect  estimate  of  e  was  made  from  the 
friction  velocity  and  the  logarithmic  wind  profile 
(Figure  6) .   Thus , 


Z  +  Z 

U*  =  KU/£N  — = 9-  (17) 

^0 


where  K  =  0.4  is  the  von  Karman  constant  and  Z~  the  rough- 
ness length.   The  relation 


e    =  U*3/KZ  (18) 


was  then  used  to  evaluate  e . 

b.  Inertial  Range  Method 

Once  the  spectrum  <$>(£)    was  determined,  e  was 
estimated  by  choosing  values  of  <J>(f)  from  the  best  straight 
line  in  the  frequency  region  in  which  a  -5/3  form  was 
observed  and  then  substituting  into  Eq.  13.   A  value  of 
0.5  was  assumed  for  the  constant  K'  based  on  results 
presented  by  Boston  (1970)  and  Williams  (1974)  . 

c.  Direct  Method 

The  kinetic  energy  dissipation  rate  was 
estimated  directly  by  summing  of  the  frequency  range 
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Figure  6.   Wind  Profiles  Runs  5  and  6 
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from  10Hz  to  2KHz  the  scales  contributing  to  viscous  dissi- 
pation according  to  Eq.  14. 

Table  II  summarizes  the  results  obtained 
from  these  three  methods.   All  values  of  e  determined  by 
the  three  methods  were  felt  to  be  representative  of  the 
actual  value.   Therefore  an  average  value  for  each  run 
at  a  given  height  was  utilized  in  the  further  analysis 
of  the  temperature  spectrum. 

2 .   Normalization  of  Spectral  Results 

a.  Velocity  Spectra 

The  velocity  spectra  were  normalized  according 
to  Eq.  3  so  that  a  comparison  could  be  made  between  dif- 
ferent runs  and  previous  experiments.   Figures  7  and  8 
are  normalized  velocity  spectra  from  Runs  5  and  6  at  the 
2m  height. 

These  figures  show  that  the  spectra  nearly 
overlap  with  a  characteristic  break  from  the  -5/3  region 
near  log  k/k  =  -1. 

b.  Energy  Dissipation  Spectra 

The  energy  dissipation  spectra  were  normalized 
according  to  Eq.  4.   A  linear  plot  was  used  to  display 
results  (Figure  9) .   Maximum  dissipation  occurred  near 
the  value  of  k/k  =  0.1  which  is  in  agreement  with  previous 
research  [Boston,  1970]  . 
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Figure  7.   Normalized  NPS  Velocity  Spectrum  Run  5  at  2m 
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Figure  8.   Normalized  NPS  Velocity  Spectrum  Run  6  at  2 
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Figure  9.   Normalized  Energy  Dissipation  Spectrum 
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B.   TEMPERATURE  SPECTRUM  RESULTS 
1.   Calculations  of  eQ 


efl  values  were  estimated  from  Eq .  16  and  the 
results  are  included  in  Table  III. 
2.   Calculations  of  K  ' 


The  values  of  K  '  were  determined  from  Eq.  15. 
Values  of  $fi(f)  were  chosen  from  the  best  straight  line 
in  the  frequency  region  in  which  a  -5/3  form  was  observed. 
The  value  of  e  used  was  the  average  value  for  a  given  run 
at  a  particular  height.   Since  e  enters  into  Eq..  15  to  the 
-1/3  power,  errors  in  estimating  z   do  not  seriously  affect 

The  results  of  the  KQ  *  calculations  are  shown  in 
Table  III.   An  average  value  of  0.88  was  obtained  which  is 
in  quite  good  agreement  with  0.81  by  Boston  (1970).   It 
was  noted  that  although  no  inherent  differences  in  tempera- 
ture spectral  shape  were  observed  between  OSU  and  NPS 
sensor  results,  the  values  of  K  '  determined  at  the  23m 
height  from  the  OSU  sensors  were  very  near  the  values 
reported  by  Williams  (1974)  at  OSU. 

3 .   Normalization  of  Spectral  Results 

a.   Temperature  Spectra 

The  temperature  spectra  were  normalized 
according  to  Eq.  8.   The  temperature  spectra  from  each 
run  at  each  height  were  normalized  and  plotted  (Figures 
10-15).   A  composite  spectrum  was  obtained  by  superimposing 
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Figure  10.   Normalized  NPS  Temperature  Spectrum  Run  5  at  2 
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Figure  11.   Normalized  NPS  Temperature  Spectrum  Run  6  at  2it 
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Figure  12.   Normalized  NPS  Temperature  Spectrum  Run  5  at  7 
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Figure  14.   Normalized  OSU  Temperature  Spectrum 

Run  5  at  23m 
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all  the  spectra  on  a  single  plot  (Figure  16) .   The  break 
from  the  -5/3  region  occurs  just  beyond  log  k/k  =  -1. 
These  spectra  show  the  shape  of  the  one-dimensional 
temperature  spectra  in  air  beyond  the  -5/3  region  v/hich 
was  one  of  the  primary  objectives  of  this  thesis.   The 
slope  of  the  spectrum  beyond  the  -5/3  break  tended  to  be 
-32/5  which  agrees  well  with  Heisenberg  (1948)  who  demon- 
strated that  a  slope  of  -7  should  exist  for  very  large 
values  of  k.   With  the  exception  of  the  temperature 
spectra  at  2m,  there  is  very  little  scatter  in  the  data 
at  high  wave  numbers . 

Figure  17  is  a  composite  plot  of  characteristic 
temperature  and  velocity  spectra.   It  is  noted  that  there 
is  very  little  difference  between  the  two  spectra  but  the 
temperature  spectrum  extends  to  slightly  higher  wave 
numbers  than  the  velocity  spectrum.   A  -32/5  slope  is 
observed  for  both  temperature  and  velocity, 
b.   Temperature  Dissipation  Spectra 

The  temperature  dissipation  spectra  were 
normalized  according  to  Eq.  9  and  again  plotted  linearly. 
Figure  18  shows  a  typical  normalized  temperature  dissipa- 
tion spectrum.   There  is  very  little  difference  between 
the  temperature  and  velocity  dissipation  spectra.   The 
normalized  shapes  are  almost  identical  with  the  temperature 
spectrum  dropping  off  slightly  faster  than  the  velocity 
spectrum. 
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Figure  16.   Composite  Normalized  Temperature  Spectra 
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Figure  18.   Normalized  Temperature  Dissipation  Spectrum 
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VI.   SUMMARY  AND  CONCLUSIONS 

The  two  primary  objectives  of  this  research  were  to 
determine  the  shape  of  the  high  wave  number  one-dimensional 
temperature  spectrum  and  the  evaluation  of  the  Kolmogorov 
scalar  constant  K  ' . 

As  seen  in  the  results  section,  this  investigator  has 
determined  that  the  high  wave  number  one-dimensional 
temperature  spectrum  falls  off  with  a  slope  of  -32/5  beyond 
the  -5/3  break  which  occurs  slightly  beyond  log  k/k   =  -1. 
Also  the  determination  of  the  Kolmogorov  scalar  constant, 
Kfi ' ,  to  be  0.88  is  considered  to  be  significant  not  only 
in  absolute  value  but  also  that  it  increases  the  data 
base  of  high  frequency  temperature  fluctuation  analysis. 
Further,  the  criticism  that  previous  evaluations  were  based 
on  records  that  were  too  short  is  removed.   It  was  also 
observed  (Table  III)  that. K  '  increased  in  value  both  with 
height  and  Reynolds  Number. 

As  an  additional  result  of  this  research,  it  is  con- 
cluded that  the  spectra  of  temperature  and  velocity  at 
high  wave  numbers  are  very  similar  with  the  temperature 
spectra  extending  to  slightly  higher  wave  numbers  than 
the  velocity  spectra  but  that  they  both  fall  off  at  the 
same  rate. 


52 


BIBLIOGRAPHY 


Batchelor,  G.  K.,  "Small  Scale  Variations  of  Convected 

Quantities  Like  Temperature  in  a  Turbulent  Fluid  (Part  1)." 
J.  Fluid  Mech.  5,  p.  113-133,  1959. 

Batchelor,  G.  K.,  Howells,  J.  D.,  and  Townsend,  A.  A., 
"Small  Scale  Variations  of  Convected  Quantities  Like 
Temperature  in  a  Turbulent  Fluid  (Part  2)."   J.  Fluid 
Mech.  5,  p.  134-139,  1959. 

Boston,  N.  E.  J.,  An  Investigation  of  High  Wave  Number 

Temperature  and  Velocity  Spectra  in  Air,  Ph.D.  Thesis , 
University  of  British  Columbia,  1970 . 

Corrsin,  S.,  "On  the  Spectrum  of  Isotropic  Temperature 

Fluctuations  in  an  Isotropic  Turbulence."  J .  Appl .  Phys . 
22,  p.  469,  1951. 

Naval  Postgraduate  School  Technical  Report  NPS-588Bb 

72021,  A  High  Frequency  Platinum  Thermometer  System 
for  Measuring  Turbulent  Atmospheric  Temperature 
Fluctuations,  by  N.  W.    J~.    Boston  and  E.  IT.    Sipe,  19  75. 

Gibson,  C.  H. ,  and  Schwarz,  H.  W. ,  "The  Universal  Equilibrium 
Spectra  of  Turbulent  Velocity  and  Scalar  Fields." 
J.  Fluid  Mech.  16,  p.  365-384,  1963. 

Heisenberg,  W.  Z.,  Physik,  124,  628,  1948. 

Lumley,  J.  L.,  and  Panofsky,  H.  A.,  The  Structure  of 
Atmospheric  Turbulence,  Wiley,  1964 . 

Nye,  J.  0.,  and  Brodkey,  R.  S.,  "The  Scalar  Spectrum  in 
the  Viscous-Convective  Subrange."  J.  Fluid  Mech.  29, 
p.  151-163,  1967. 

Panofsky,  H.  A.,  "The  Spectrum  of  Temperature."  Radio 
Science,  4,  p.  1143-1146,  1969. 

Stewart,  R.  W. ,  Wilson,  J.  R. ,  and  Burling,  R.  W. ,  "Some 
Statistical  Properties  of  Small  Scale  Turbulence  in  an 
Atmospheric  Boundary  Layer."  J.  Fluid  Mech.  41,  Part  1. 
p.  141-152,  1970. 

Williams,  R.  W.,  Jr.,  High  Frequency  Temperature  and 
Velocity  Fluctuations  in  the  Atmospheric  Boundary 
Layer,  Ph.D.  Thesis ,  Oregon  State  University,  Corvallis, 
1974. 


53 


Wyngaard,  J.  C,  "The  Effect  of  Velocity  Sensitivity  on 
Temperature  Derivative  Statistics  in  Isotropic 
Turbulence."  J.  Fluid  Mech.  48,  p.  763-769,  1971. 


54 


INITIAL  DISTRIBUTION  LIST 

No.  Copies 

1.  Defense  Documentation  Center  2 
Cameron  Station 

Alexandria,  Virginia    22314 

2.  Library,  Code  0212  2 
Naval  Postgraduate  School 

Monterey,  California   93940 

3.  Department  Chairman,  Code  58  3 
Department  of  Oceanography 

Naval  Postgraduate  School 
Monterey,  California   93940 

4.  Assoc  Professor  N.  E.  J.  Boston,  Code  58  Bb    1 
Department  of  Oceanography 

Naval  Postgraduate  School 
Monterey,  California   93940 

5.  Assoc  Professor  T.  M.  Houlihan,  Code  59  Hm    1 
Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California   93940 

6.  Commanding  Officer  1 
Fleet  Numerical  Weather  Central 

Monterey,  California   93940 

7.  Commanding  Officer  1 
Environmental  Prediction  Research  Facility 

■  Monterey,  California   93940 

8.  Department  of  the  Navy  1 
Commander  Oceanographic  System,  Pacific 

Box  139  0 

FPO  San  Francisco   96610 

9.  Oceanographer  of  the  Navy  1 
Hoffman  II 

200  Stovall  Street 
Alexandria,  Virginia   22332 

10.   Office  of  Naval  Research  1 

Code  480 
Arlington,  Virginia   22217 


55 


11.  Dr.  Robert  E.  Stevenson 
Scientific  Liaison  Office,  ONR 
Scripps  Institution  of  Oceanography 
La  Jolla,  California   92037 

12.  Library,  Code  3330 

Naval  Oceanographic  Office 
Washington,  D.  C.   20373 

13.  SIO  Library 

University  of  California,  San  Diego 

P.  0.  Box  2367 

La  Jolla,  California   92037 

14.  Department  of  Oceanography  Library 
University  of  Washington 
Seattle,  Washington    98105 

15.  Department  of  Oceanography  Library 
Oregon  State  University 
Corvallis,  Oregon   97331 

16.  Dr.  Carl  Gibson 

University  of  California,  San  Diego 

Department  of  Ames 

P.  0.  Box  109 

La  Jolla,  California   92037 

17.  Dr.  James  J.  O'Brien 
Program  Director 
Physical  Oceanography 

Ocean  Science  and  Technology  Division 
Office  of  Naval  Research 
Arlington,  Virginia   22217 

18.  Dr.  R.  M.  Williams 
Oregon  State  University 
School  of  Oceanography 
Corvallis,  Oregon   97331 

19.  Dr.  J.  Wyngaard 
AFCRL  (LYB) 

L.  G.  Hanscom  Field 

Bedford,  Massachusetts   01730 

20.  LT  Robert  Thomas  Simril,  USN 
2321  Wensley  Drive 
Charlotte,  North  Carolina   28210 


56 


21.  Dr.  R.  w.  Burling 
Institute  of  Oceanography 
University  of  British  Columbia 
Vancouver  8,  British  Columbia 
Canada 

22.  Dr.  Michel  Coantic 

Institut  de  Mecanique  Statistique 
de  la  Turbulence 
12,  Avenue  General  Leclerc 
Marseille  (3e)  ,  France 

23.  Dr.  A.  Gyr 

Institute  of  Hydromechanics  and 
Water  Resources  Management 
8006  Zurich,  Tannens trasse  1 
Switzerland 

24.  Dr.  Sjrfren  Larson 
Research  Establishment  Risrf 
DK-4000  Roskilde,  Denmark 


57 


Thesis 

j  dL 

S517 

S  imri  1 

c.l 

High   frequency   tem- 

perature  fluctuations 

in    the   atmospheric 

boundary    layer. 

Thesis  .32 

S517    S  imri 1 

c.l  High   frequency   tem- 

pera tv  re    f  l"ct"ations 
in   the  atmospheric 
boundary    layer. 


thesS517 

H!^lTle™ltemP™™  fluctuations 

fl'f'l'lll1'!!' I'lll  Illl  Illl 


-««ll  Hill  lllll  ill 

3  2768  001  91435  1 

DUDLEY  KNOX  LIBRARY 


r 


I 


^m 


■ 

■ 

■Mid 

m 


